Fluctuating initial condition and smoothening effect on elliptic and triangular flow 
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In heavy ion collisions, event-by-event fluctuations in participating nucleon positions can lead to 
triangular flow. Generally, one uses Monte-Carlo Glauber model to obtain the participating nucleon 
positions. To use in a hydrodynamic model, the positions needs to be smoothened. We study the 
effect of smoothening of Glauber Monte-Carlo initial conditions on elliptic and triangular flow. It 
is shown that integrated as well as differential elliptic and triangular flow remain largely unaltered, 
irrespective of functional form of the smoothening function, or the smoothening parameter 
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Azimuthal distribution of the produced particles is one 
of the important observables in ultra-relativistic nuclear 
collisions. In a non-zero impact parameter collision be- 
tween two identical nuclei, the collision zone is asym- 
metric. Multiple collisions transform the initial asymme- 
try into momentum anisotropy. Momentum anisotropy 
is best studied by decomposing it in a Fourier series, 



dN 



N 
27 



2 v n cos(n(f> 



nipn) 



1,2,3. 



(1) 

(j> is the azimuthal angle of the detected particle and ip n 
is the plane of the symmetry of initial collision zone. 
For smooth initial matter distribution, plane of sym- 
metry of the collision zone coincides with the reaction 
plane (the plane containing the impact parameter and 
the beam axis), ip n = ^/jp,Vn. The odd Fourier co- 
efficients are zero by symmetry. However, fluctuations 
in the positions of the participating nucleons can lead 
to non-smooth density distribution, which will fluctuate 
on event-by-event basis. The participating nucleons then 
determine the symmetry plane (tppp), which fluctuate 
around the reaction plane As a result odd harmon- 
ics, which were exactly zero for smoothed initial distri- 
bution, can be developed. It has been conjectured that 
third hadronic U3, which is response of the initial tri- 
angularity of the medium, is responsible for the observed 
structures in two particle correlation in Au+Au collisions 
0,0,[3,[!,[1,0- The ridge structure in pp collisions 
also has a natural explanation if odd harmonic flow de- 
velops. Recently, ALICE collaboration has observed odd 
harmonic flows in Pb+Pb collisions Q. In most central 
collisions, the elliptic flow (^2) and triangular flow (1)3) 
are of similar magnitude. In peripheral collisions how- 
ever, elliptic flow dominates. 

In a hydrodynamic model, collective flow is a response 
of the spatial asymmetry of the initial state. For ex- 



ample, elliptic flow is the response of ellipticity of the 
initial medium, triangular flow is the response of the ini- 
tial triangularity of the medium and so on. In theoreti- 
cal simulations, one generally uses Monte-Carlo Glauber 
model to obtain the event-by-event initial conditions. In 
a Monte-Carlo Glauber model, according to the density 
distribution of the colliding nuclei, two nucleons are ran- 
domly chosen. If the transverse separation between them 
is below a certain distance, they are assumed to inter- 
act. Transverse position of the participating nucleons is 
then known in each event and will fluctuate from event- 
to-event. If a particular event has N part participants, 
participants positions in the transverse plane can be la- 
beled as, (xi,yi),(x2,y2)-—(x n part,Vnpart)- The energy 
density in the transverse plane can be approximated as, 
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However, fluid dynamical model require continuous 
density distribution and discrete distribution as in EqJ5] 
cannot be evolved. To use in a hydrodynamic model, the 
discrete density distribution has to be converted into a 
smooth energy-density distribution. This can be done by 
smearing the discrete participating nucleon position by 
some smoothening function, S(x — Xi,y — yi) — > g{x — 
Xi, y — yi, Ci, & •■)■> d being parameters of the smoothen- 
ing function g. Incidentally, even though smoothening of 
the nucleon positions is a must if one uses Monte-Carlo- 
Glauber model, existing literature in general do not elab- 
orate on the procedure. However, if the flow coefficients 
are to be used for diagnostic purpose, e.g. to constrain 
the viscosity over entropy ratio, it is important to know 
the effect of smoothening. In the present brief report, we 
study the effect of smoothening of Monte-Carlo Glauber 
model initial condition on the elliptic and triangular flow. 
It will be shown that the flow coefficients depend mini- 
mally on the smoothening function. 

In the following, we write the smoothen out energy 
density in the transverse plane as, 
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where g{x — Xi, y — yi, Ci, C2—) is the smoothening func- 
tion and Q's are the parameters of the smoothening func- 
tion. So is a parameter, which is fixed such that event 
averaged particle multiplicity reproduces the experimen- 
tal value. In the following, we consider two smoothen- 
ing functions, (i) a Woods-Saxon distribution and (ii) a 
Gaussian distribution. The Woods-Saxon distribution is, 
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The half radius C in EqH]is kept fixed C — lfm. Only 
the diffuseness parameter a is varied to obtain different 
smoothening. In the following, we consider four differ- 
ent values, a=0.05, 0.1, 0.25 and 0.5 fm. The effect of 
smoothening on the energy density distribution can be 
seen in FigQ] In panel (a) , the transverse distribution of 
participating nucleons, in b=6.5 fm Au+Au collisions, for 
a typical Monte-Carlo event is shown. The discrete posi- 
tions are smoothened with a Woods-Saxon distribution. 
In panel (b)-(d), the smoothened out energy density dis- 
tribution, for diffuseness parameter a=0.1, 0.25 and 0.5 
fm are shown. One notes that structures in the density 
profiles are diffused as the diffuseness increase. In the 
following, we study the effect of smoothening on elliptic 
and triangular flow. 

To show the dependence on the functional form, we 
have also smoothened the participant nucleons positions 
with a Gaussian function, 
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Two values of the Gaussian width, cr=0.1 and 0.5 fm 
are considered. In the following, we will refer to the dif- 
fuseness parameter (a)of the Woods-Saxon function and 
the width (a) of the Gaussian as the smoothening pa- 
rameter (s). 

With fluctuating initial conditions, in each event, 
space-time evolution of the fluid is obtained by solving 
the energy-momentum conservation equations d^T^" — 
0, in (t = yjt 2 — z 2 , x, y, ij s — ^lnf^f) coordinate sys- 
tem, with the code AZHYDRO-KOLKATA 0. We as- 
sume boost-invariance. We disregard any dissipative ef- 
fect. Ideal hydrodynamics equations are closed with an 
equation of state (EoS) p = p(e). Currently, there is con- 
sensus that the confinemcnt-dcconfinement transition is 
a cross over and the cross over or the pseudo critical tem- 
perature for the transition is T c «170 MeV [j"oWl3l|. In 
the present study, we use an equation of state where the 
Wuppertal-Budapest [l(| [H| lattice simulations for the 
deconfined phase is smoothly joined at T = T c = 174 
MeV, with hadronic resonance gas EoS comprising all 
the resonances below mass in res =2.5 GeV. Details of the 
EoS can be found in [bH ]. 

In addition to the initial energy density for which we 
use the smoothened out Monte-Carlo Glauber model, 
solution of hydrodynamic equations requires to specify 
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FIG. 1: (color online) (a) discrete positions of participat- 
ing nucleons in the transverse plane in a typical event from 
Monte-Carlo Glauber model, in b=7 fm Au+Au collision, (b)- 
(d) density distribution obtained by smoothening the discrete 
position by Woods-Saxon function, with diffuseness parame- 
ter, o=0.1, 0.25 and 0.5 fm. C = lfm. 



the thermalisation or the initial time t; and fluid veloc- 
ity (v x (x,y),v y (x,y)). A freeze-out prescription is also 
needed to convert the information about fluid energy den- 
sity and velocity to invariant particle distribution. We as- 
sume that the fluid is thermalized at Tj= 0.6 fm and the 
initial fluid velocity is zero, v x {x,y) — v y (x,y) = 0. The 
frcczc-out is fixed at Tp—130 MeV. We use Cooper-Frye 
formalism to obtain the invariant particle distribution 
from the freeze-out surface. From the invariant distribu- 
tion, harmonic flow coefficients are obtained as [15] . 
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In a boost-invariant version of hydrodynamics, flow co- 
efficients are rapidity independent and in the following. 
Present simulations are suitable only for central rapidity, 
y «0, where boost-invariance is most justified. Hereafter, 
we drop the rapidity dependence. tp n in Eg 16I7I is the re- 
action plane angle. We use the following definitions of 
initial eccentricity £2 and triangularity £3 ;S] , 6] , , 
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In the following we show simulation results for 30- 
35% Au+Au collisions at ^/s NN —200 GeV. We have con- 
sidered N event = 100 events. Recently in fl6| . it was 
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FIG. 3: (color online) Top panel: integrated elliptic flow 
scaled by initial eccentricity as a function of smoothening 
parameter (s). The filled and solid symbols are for Woods- 
Saxon and Gaussian smoothening function respectively. Bot- 
tom panel: same for the initial triangularity scaled integrated 
triangular flow. 



FIG. 2: (color online) Charged particles multplicity (top 
panel), integrated elliptic flow (the middle panel) and in- 
tegrated triangular flow (the bottom panel), as function of 
smoothening parameter (s). The filled and solid symbols are 
for Woods-Saxon and Gaussian smoothening function respec- 
tively. 

shown that with fluctuating initial conditions, event av- 
eraged as well as variance of elliptic flow and triangu- 
lar flow remain approximately unaltered for N event =50- 
2500. N event = 100 is then sufficiently large to study 
the effect of smoothening on flow coefficients. The pa- 
rameter £rj in EqJ21 is fixed such that irrespective of the 
smoothening function or smoothening parameter, hydro- 
dynamic evolution approximately reproduces the experi- 
mental charged particles multiplicity. Results of our sim- 
ulations for N — 100 Monte-Carlo Glauber model events 
are shown in FigJ^El The filled/open symbols are event 
averaged values and the bars represent the variance. The 
filled symbols are obtained with smoothening with the 
Woods-Saxon function, the open symbols are obtained 
with the Gaussian function. The top panel of FigJ2]shows 
that the simulations are constrained such that the aver- 
age charged particles multiplicity is dN/dy 210, irre- 
spective of the smoothening function and the smoothen- 
ing parameter. In the middle panel of FigJ^J depen- 
dence of the integrated elliptic flow on the smoothening 
function and smoothening parameter is shown. Both for 
Gaussian and Woods-Saxon smoothening function, ellip- 
tic flow do not show any dependence on the smoothening 
parameter. Note that for the Woods-Saxon function, the 
smoothening parameter is varied by a factor of 10 and for 



the Gaussian function, by a factor of 5. Even though the 
smoothening parameter is changed a factor of 5-10, event 
averaged vi and their variances remain approximately 
unchanged. More interestingly, smoothening with Gauss 
or Woods-Saxon function, both produces approximately 
same i>2- 

In the botom panel of FigJ2J simulation results for the 
triangular flow is shown. As it was observed for the el- 
liptic flow, integrated triangular flow also do not show 
any dependence on the functional form of the smooth- 
ing function or on the smoothening parameter. Both for 
Gaussian function and Woods-Saxon function, approxi- 
mately similar value is obtained for the triangular flow. 
Event averaged triangular flow also show marginal depen- 
dence on the smoothening parameter. In event-by-event 
hydrodynamics, with Monte-Carlo Glauber model initial 
conditions, the elliptic and triangular flow seems to be 
independent of the smoothening function as well as of 
the smoothening parameter. 

In ideal hydrodynamics, anisotropic flow is expected to 
scale with initial asymmetry of the reaction, quantified in 
terms of initial eccentricity and initial triangularity etc. 
Presently, we donot study the scaling behavior of flow co- 
efficient. Rather, we study the effect of smoothening on 
the scaled observables, the initial eccentricity scaled ellip- 
tic flow [vijt-i) and initial triangularity scaled triangular 
flow («2/es)< The simulation results are shown in Fig|3] 
As it was for the elliptic and triangular flow, initial eccen- 
tricity scaled elliptic flow as well as initial triangularity 
scaled triangular flow also remain unaltered for Gaussian 
and Woods-Saxon smoothening functions. It also show 



4 



0.5 
0.4 
,0.3 
0.2 
0.1 

0.5 
0.4 
0.3 
0.2 
0.1 




v s = .05 fm 
■ □ s = .10 fm 

a s = .25 fm 
+ * s = .50 fm 



if 



0.5 



1 1.5 2 

P T (GeV/c) 



2.5 



FIG. 4: (color online) Dependence of differential elliptic (the 
top panel) and differential triangular (the bottom panel) 
flow, on the smoothening parameter (s). The open and 
filled symbols are respectively for Woods-Saxon and Gaus- 
sian smoothening function. 



little dependence on the smoothening parameter. 

Lastly, in FigJH effect of smoothening on the differ- 
ential elliptic and triangular flow is shown. As before, 
the filled/open symbols are for Woods-Saxon/Gaussian 
smoothening functions. As it was for the integrated V2, 
differential vi also depend marginally on the smoothening 
function or on the smoothening parameter. For Gaussian 
smoothening function, event averaged v^lpr ~ 1.5GeV) 



changes by ~ 3% for a factor of 5 change in the smoothen- 
ing parameter, from 0.1 to 0.5. The change is even less 
for the Woods-Saxon function. Indeed, within the un- 
certainties, the differential vi is approximately indepen- 
dent of the smoothening function and of smoothening pa- 
rameter. The differential triangular flow is also remains 
approximately independent of the smoothening function 
and smoothening parameter. However, absolute change 
is marginally more. For example, for Gaussian smoothen- 
ing function, event averaged V2(j>t ~ 1.5GeV) is reduced 
by ~ 10% for a factor of 5 increase in the smoothening 
parameter. The change is even less (~7%) for a factor of 
10 change in the smoothening parameter of the Woods- 
Saxon function. 

To summarise, in event- by- event hydrodynamics, 
discrete participant nucleon position from Monte- Carlo 
Glauber model calculation required to be smoothened. 
We have investigated the effect smoothening on elliptic 
and triangular flow. In each event, participating nucleons 
position from Monte-Carlo Glauber model calculations 
are smoothened with either a Gaussian function or a 
Woods-Saxon function. Smoothening is controlled by 
a smoothening parameter, e.g. diffuseness for Woods- 
Saxon function and width for the Gaussian function. It 
was shown that integrated as well as differential elliptic 
and triangular flow remain largely unaltered, irrespective 
of functional form of the smoothening function, or the 
smoothening parameter. 
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